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Abstract 

The Hamiltonian formalism of bigravity and bimetric theories is 
constructed for the general form of potential between two metrics. It 
is natural to study the role played by the lapse and shift functions in 
theories with two metrics on the base of Kuchar's approach because 
then they do not depend on the choice of space-time coordinate frame. 
The algebra of first class constraints is derived in Dirac brackets con- 
structed from the second class constraints. It is the celebrated algebra 
of the hypersurface deformations. Fixing one of the metrics we ob- 
tain a bimetric theory without any first class constrains. Then we can 
use the symmetries of the background metric to construct conserved 
quantities looking ultralocally when written through the potential. 
The special case of potential providing the less number for degrees of 
freedom will be treated in the companion paper. 

1 Introduction 



The idea to use more than one dynamical metric for the description of the 
real Universe is not rather new. This theory is usually called multi-gravity. 
There is a hope that such a modification of the gravitational theory may allow 
to solve the dark energy problem. We restrict ourselves to the two metrics 
case, and will call such a model as bigravity, following pioneers pp. Suppose 
each metric interacts only with its own sort of matter, but these metrics can 
interact with each other by means of some potential, i.e. a scalar density 
algebraically constructed of the two metrics. If we fix one metric by giving 
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it the absolute meaning, we move from bigravity to bimetric theory, as since 
N. Rosen's [2] times the metric theories of gravity on the fixed background 
are called. In particular, massive gravity theories [3] applied for explaining 
the dark energy phenomena in the Universe are in this category. 

We plan to consider the Hamiltonian formalism of bigravity in two papers 
dividing our presentation in two parts which we call the general case and the 
special case. In general case the concrete form of interaction potential is 
not given, we only say that it is constructed algebraically (ultralocally, i.e. 
without derivatives), it is a scalar density and does not contain any fields of 
matter. We call the special case the model with a potential giving one degree 
of freedom less than the general case. The form of such a potential had been 
proposed recently [I]. Unfortunately straightforward calculations with this 
potential in Hamiltonian formalism are difficult because of its matrix square 
root form, and so the analysis given already in a publications [5] does not 
look transparent and decisive. We are intended to derive the requirements 
on the potential which are necessary for realizing the program proclaimed in 
papers [3], starting from a potential of the general form. In contrast to the 
previous articles on this subject we exploit not ADM, but Kuchar's approach 
which preserves the freedom of spacetime coordinates choice. 

We denote as f^ v and call as the first metric the metric that will be 
fixed (i.e. become a background one) after a shift from bigravity to bimetric 
theory, and the second metric g^ u which will always stay dynamical. Greek 
indices will run from to 3, Latin - from 1 to 3, signature is (— , +, +, +). 

Consider two copies of General Relativity Lagrangian, each copy supplied 
with its own sort of matter as a source (interaction with matter is supposed 
minimal) , 



where f^ u and g^ v are the first and the second spacetime metrics correspond- 
ingly, / and g are their determinants, Rv) and Rffi are Ricci tensors, 

and are gravitational constants, dfj and are Lagrangians of the 
first and the second matter, A denotes abstract indicies for matter fields, 
and construct of them a new Lagrangian by adding a potential 




(1) 






(3) 
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therefore, 



C= CP + + J~fU(f^,g^). 



(4) 



It is clear that a choice of metric determinant which is present explicitly is 
not of importance because function U(f lu/ ,g fiV ) may arbitrarily depend on 
the ratio of two determinants. 

2 Kuchaf's Hamiltonian formalism in Gen- 
eral Relativity 

Let the action of General Relativity (GR) is written as follows 



where X a are spacetime coordinates, and the form of Lagrangian is the same 
as in formulas (DO), (El). 

In the construction of Hamiltonian formalism it is necessary to separate 
explicitly the time coordinate from the spatial ones. The state should be de- 
termined by values of the gravitational and matter fields given at all points 
of space at one definite moment of time. So, the state is to be given on a 
spacelike hypersurface embedded in spacetime. The flow of time corresponds 
to the motion of this hypersurface through spacetime or to the continuous 
transform from one hypersurface to another, so we need one-parametrical 
family of spacelike hypersurfaces. The role of time coordinate t is to be 
played by a parameter which continuously and monotonically numerates hy- 
persurfaces. It is suitable to introduce spatial coordinates x l on the initial 
hypersurface and then continue them to all hypersurfaces in such a way that 
lines going through the points with the same values of coordinates can be 
treated as observer world lines, i.e. they should be timelike. 

In ADM approach [6], which is the most popular, the choice of spacelike 
hypersurfaces family is determined by the choice of spacetime coordinate 
frame: 



t = X°, x* = X\ 7y = (fc, N=(-g 00 )- 1 / 2 , Ni = goi . (6) 



There is another approach proposed by Kuchar [7] where two coordinate 
systems are in action, one is an arbitrary spacetime coordinate system X a , 




(5) 



3 



the other (t, x l ) is defined by embedding variables 

X a = e a (x\t), (7) 

where time is a parameter monotonically numerating hypersurfaces and other 
three coordinates numerate points on a hypersurface. We will follow this fully 
covariant Kuchar's method. Then fields 

de a 

will simultaneously be vectors in spacetime and covectors in space. The 
metric induced on a hypersurface is given as follows 

lij 9fJ,v&i 6 j ■ (9) 

As usual, inverse matrices to metrics g^ u and 7^ are denoted as g^ v and ^ . 
Moving indices (both Greek and Latin) up and down is provided by means 
of them: 

e ai = gape?, e l a = g aP e^ K '. (10) 

We will use here a bar to make a distinction between the quantities defined by 
means of metric and the quantities ef which are independent of the metric. 
In the following, when we will work with two metrics, the bar will stay 
connected with g^ u . 

Next we introduce normal 1-form 

n a e? = 0, (11) 

which can be converted with the help of metric tensor to a normalized vector: 

n a =g a %, = -1. (12) 

We can use basis (n a , ef) to decompose every vector or tensor in spacetime, 
for example, Ricci tensor (with account for its symmetric nature): 

W = RU-rfn" + R lL {e ! tn u + n^) + We^e), (13) 

here the components are calculated as follows: 

R ±± = R^n^n u , R l± = -ETn^e*, R ij = R^e^e vj . (14) 
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With the known technique [7J [S], based on expression of the Riemann 
tensor projections as the corresponding projections of covariant derivatives 
commutators, it is possible to write GR Lagrangian density containing the 
single spacetime metric g^ v in the following form 

a = ^ (J*M - K 2 + SpK 2 ) + £ M (0 A , N, N*) t (15) 

where boundary terms, i.e. total time derivatives and spatial divergences are 
ignored, as we do not discuss boundary conditions. Here = 1/16ttG^ 9 \ 
is the scalar curvature of the induced metric, 7 = det||7y||, lapse and 
shift functions N, N % are components of the decomposition of time vector 
field over the basis constructed with metric g^ v 

BX a 

N a = — — = Nn a + N*ef. (16) 

The second fundamental form Kij arising in the Lagrangian density can be 
expressed through the already introduced variables by formula: 

^; = ^(% + %-7y)> (17) 

here K = ^K^, Sp^ 2 = K^K 1 ^ , vertical bar denotes a covariant derivative 
defined by induced metric. The momenta conjugate to variables 7^ and <p A 
are determined as follows: 

and the momenta conjugate to N and iV* are zero, as expression (ITB1) does 
not contain velocities N, N 

7Tjv = 0, 7Tjvi = 0, (19) 

so, these equations are primary constraints. The velocities for gravitational 
variables 7^ can be expressed by means of momenta 

la = + N j{i + —j=-faj ~ 7,, 2 )• (2°) 
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We believe that it is possible to do the same for the matter variables (in other 
case we can analyse new constraints arising there), and after performing the 
Legendre transform 

H canonica i = J d 3 x (jr^iij + tt a< P a - C) , (21) 

to obtain the canonical Hamiltonian (up to surface terms) 

H canonical = J d 3 x (NH + N l n t ) . (22) 

Poisson brackets for functionals of canonical variables 21X6 cLS follows 

rr/11 f j3 (5F5G 6F 5G 5F 5G 5F 5G /rn ^\ 

{F, G} = d 3 x[ — — + — + ^^ + tttt f H G . 

1 J J \5N5tt n 5N 1 5tt Ri &y ij &it % 3 5<f) A Sir A v ') 

(23) 

According to Dirac method [9] , one should further go to the extended Hamil- 
tonian by adding the primary constraints with arbitrary Lagrangian multi- 
pliers 

Hextended = / <Px (NH + + \tT N + AV^) . (24) 



Then expressions 



«^^(^ (7>+K<J, (y- Spi2 )) +flM 



(25) 



and 

H i = -2ni\ j + 'HiM (26) 

become secondary constraints, as their equality to zero is necessary for pri- 
mary constraints staying zero in evolution. Here 7r = jijir^, Spn 2 = TT^Xy, 
'H-Mi'fiiM is the matter contribution. 

The standard procedure requires to remove four pairs of canonical vari- 
ables N,N\ 71^,71^ and of primary constraints (fT9j) by putting on four 
gauge conditions which really change only the meaning of four letters, i.e. 
replace canonical variables with functions (we do not like to introduce new 
notations and we make this change only here) 

N-tf(x) = 0, N i -AT i (x) = 0. (27) 
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These gauges have nonzero Poisson brackets with primary constraints (1191) 
and allow to treat all of them together as eight second class constraints, this 
leads to Dirac brackets: 

The number of gravitational degrees of freedom is now determined by a simple 
calculation: take the number of induced metric 7^ independent components 
and subtract the number of first class constraints: 6—4 = 2, now Hamiltonian 
takes a following form 

Hpartially reduced J d 3 x{Afn+Af l ni) . (29) 

Avoiding cumbersome terminology and notations we as usual will call Dirac 
brackets fl28j) as Poisson brackets, the partially reduced Hamiltonian fT29|) as 
Hamiltonian, N = Af, N l = M l as Lagrangian multipliers standing before 
constraints ( 12"5|) . f l2"B"j) . These constraints are first class as their algebra is as 
follows 

{H(x),U(y)} = (V lk (x)n k (x) + r, ik {y)U k {y))SA^ v), (30) 
{Hi{x),Uk{y)} = Hi(y)5,k{x, y) + 'H k (x)5^(x, y), (31) 
{%{x),U{y)} = n{x)8^x,y), (32) 
reflecting a freedom of hypersurface deformations in Riemannian space. 

3 Hamiltonian approach in bigravity 

Having two spacetime metrics we get two induced metrics on a spatial hy- 
persurface: 

lij ^ji Vij f iw&i ^ j 1 (33) 

and also two different unit normal vectors that will be denoted as n a n a : 

n a e° = = n a e^ f a ^n a np = -1, g aP n a np = -1, (34) 

and two bases (n a ,e°), (n a ,e°). Moving indices up and down is provided 
by the corresponding metric tensors f^ u ,Vij an d g^vilij- We will decompose 
any spacetime vectors and tensors over basis (n a ,e° l ), constructed with the 
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help of f^v which will be called the first metric. For example, we decompose 
the second metric tensor (with account for its symmetry) as follows 

gT = g^nW + g ±l (eW + n»t%) + g^efe", (35) 

whereas, 

fi* = _ n » n v + e v e v.rf\ (36) 
the components are calculated by formulas: 

g ±± = g^n,n u , g Xi = -g^el, g* = g^e^ = 7 « + 9 -^. (37) 

We write each metric contribution to bigravity Lagrangian (j3J) in the similar 
way to equation ffT5l) . The potential, in its turn, depends on the two space- 
time metrics, therefore it is necessary to brake their symmetry and take one 
of the metrics for construction of the basis in order to transform the poten- 
tial to (3 + l)-form. Previously we already have decided to call f^ u the first 
metric and to use the basis constructed with its help. Then besides lapse N 
and shift N 1 

N a = Nn a + JV*<£ = Nn a + JV*e? (38) 

and induced spatial metrics rjij, 7 y, the potential depends on other four 
components of metric g^ v decomposition over basis (n a , ef), let us accentuate 
that in contrast to N, quantities g^^ig 1 ' 1 enter the potential in a nonlinear 
way. Relations between the two bases 

g ±i 



n a = J-g^n* - ef (39) 



-g ±L 

allow us to express lapse and shift of metric g^, i.e. N, N l through N, 
and normal projections g , g . 

N - ■ 

N = , N l = N l - —-—N. (40) 

^ 9 

It is suitable for the following to introduce new variables 

1 • a u 

ul = ~^ (41) 
-g ±A - 9 
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having the simple geometric meaning: u is an inverse of a norm (calculated 
in the second metric) of vector n a , constructed as a unit normal (in the first 
metric) to the hypersurface, and u l are three projections (calculated in the 
second metric) of coordinate basis vectors onto this unit normal 

- = -T-i— u> = Jf^k=. (42, 

Then equations (l4"0j) takes a form 

N = uN, N* = N { + u l N. (43) 

We consider as dynamic variables two sets of matter fields (pAi^A, two in- 
duced metrics on the hypersurface 7^-, components of the time vector in 
the chosen basis (i.e. lapse and shift) N, N l and four additional variables 
u,u l , taken intead of the second metric projections g ±± ,g ±l . Momenta are 
defined in the usual way 

TTM d ^ ii TT d £ d £ tAA\ 

orjij fry i: j dip A d(p A 

The two procedures of defining momenta and Legendre transform are real- 
ized in parallel and independent way for both terms & and CP . Primary 
constraints arise because velocities of eight variables N, N l , u, u 1 are absent 
in the Lagrangian density, these constraints have a following form 

TT N = 0, TC N i = 0, (45) 

vr u = 0, 7T U > = 0. (46) 

Certainly, some additional primary constraints can appear due to gauge in- 
variance which matter fields may have, but we will not discuss this here. The 
canonical Hamiltonian is derived by the Legendre transform 

H canonica l = J d 3 X (ll% + TT^ + TtyA^ + - C) (47) 

and expression of velocities through the corresponding momenta, this pro- 
vides us with the two similar looking terms without common variables plus 
the potential: 

H canonica i = / d 3 x {NH + WHi + NH + N*Hi + N , (48) 
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or in other way, 

H canonica l = j d^X (N (H + uH + li'tii + y/fjU) + JV* fa + Hi)) . (49) 

The canonical Poisson brackets have a standard appearance corresponding to 
the set of conjugate variables (77^, IT-'), (7^, 7r lJ ), (ip A , Ha), (4> A , ^a), {N, ttn), 
(N l ,7r N i), (u,tt u ), (u l ,7T u i). The requirement of conservation of the primary 
constraints (14"5]) . ( )46|) in Hamiltonian evolution forces us to put their Poisson 
brackets with Hamiltonian to zero 

{tt n ,R} = 0, {vt^,H} = 0, K,H} = 0, {tv,H} = 0, (50) 

this leads to the secondary constraints: 

TZ = H + uH + trtU + y/TjU = 0, (51) 

Ki = Hi + Hi = 0, (52) 

where U = *JrjU. In the general case which is under discussion in this pa- 
per the four last constraints together with the four last primary constraints 
occurs second class pU] and can be excluded from Hamiltonian after intro- 
duction of Dirac brackets (see Appendix [B]). The special case will be treated 
in the companion article. For other secondary constraints the algebra of 
hypersurface deformations should be fulfilled: 

{TZ(x), TZ(y)} D = [r, ik (x)TZ k (x) + rj ik (y)K k (y)] S ti (x, y), (54) 

{TZi(x),n k (y)} D = Ki{y)5 >k (x, y) + TZ k (x)5 :i (x, y), (55) 
{Ki(x),K(y)} D = K(x)5 ti (x,y). (56) 

Calculating l.h.s. first as Poisson brackets, then applying second class con- 
straints (|53|) . at last comparing results with r.hs., we derive conditions on the 
potential given as function of variables 77^, 7^, u, u l (details see in Appendix 

2 Vjk ^- + 2 ljk ^- - u'^r = SlU, (57) 



' driij ()- ,j du k 
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These conditions preserve their form under a change of variables correspond- 
ing to the interchange of bases and roles of the two metrics: 

1 • v? ~ 

u — y — , u l — y , rjij <H- 7^, U -» uU. (59) 

u u 

The simplest method to derive Hamiltonian equations is to use canonical 
Hamiltonian and Poisson brackets. Then they are looking similar to GR 
Hamiltonian equations (here nabla denotes a covariant derivative defined by 
%): 

^ = { r]lJ ,Rf} = V j N i + V i N j + ^'^(U ij -r ]i ~), (60) 

7y = { 7y ,H ff } = % + % + «C fl ^(7ry-7 i ~), (61) 

(9/7 

= {n«H/}-JV— , (62) 

7r iJ = {vr^,H 9 } -iV-— , (63) 

^ = (64) 

(65) 

but one should remember that if N, N l are arbitrary Lagrangian multipliers 
(before putting on gauge conditions), u, u l after solving second class con- 
straints become functions of which in their turn are expressed by 
equations fl25|) . fl26|) through canonical variables. Certainly, these Hamilto- 
nian equations can be derived by means of Dirac brackets from Hamiltonian 
after applying second class constraints to it: 

When we are given embedding variables X a = e a (t, x l ), we can express both 
spacetime metrics in coordinates X a , if we first estimate the unit normal 
vector 

1 fde a ■ \ 

* = n (■*-»'*)• (67) 

and then apply formulas ( l35l) . ( l36l) . ( |37l) for two metrics decompositions in 
this basis. 
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4 Hamiltonian approach in bimetric theories 



Consider now the situation when only the second metric is dynamical and 
the first one is background, i.e. it is a fixed solution of GR equations. Tra- 
ditionally [2] such theories are called bimetric. Recently they attracted a lot 
of attention and some valuable reviews have appeared [TU [121 H3J EH]- 

If metric tensor and embedding variables e a (t,x k ) are given, then we 
already know functions N(t, x k ), N l (t, x k ), r]ij(t, x k ). We can take expression 
([49]) . derived above for bigravity, as a Hamiltonian, treating N,N l not as 
variables, but as given parameters, and taking into account that are 
now zero, 

H = J d 3 x (iV (uH + vfiU + y/ijU) + NWi) . (68) 

Then primary constraints are only equations ( HBl) . and they lead to secondary 
constraints (1531) . In the general case (remember, that special case will be 
treated in another paper) all these constraints are second class, and Dirac 
brackets constructed from them will coincide with Poisson brackets for the 
functionals which depend only on the independent variables 7y, n 1 ^ (f) A , tia- 

For example, Hamiltonian fits in this category after excluding variables u, ir u , 
u\ n u i which can be done with the hekp of equations (j4~6]) . ([53j) . Just opposite, 
we can express Hamiltonian by means of the potential U and its derivatives 
with respect to variables u,u l : 

»- !*(»(*- 4 m 

in this case more complicated formulas for Dirac brackets which take into 
account variables u,u l are to be used, see Appendix [El 

Hamiltonian of bimetric theory is not a linear combination of first class 
constraints, in contrast to bigravity and GR Hamiltonians, so we can use 
it for construction of conserved quantities. For example, if the background 
metric is flat, then the Hamiltonian density may be interpreted as a density 
of energy, momentum or angular momentum, and all these quantities are 
ultralocal functions of the corresponding lapse and shift N,N\ two induced 
metrics %-,7ij and u,u l . 
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5 Conclusion 



In this article we have studied the theory of bigravity with a potential of 
the general form and built the Hamiltonian formalism, avoiding unnecessary 
noncovariance in defining the lapse and shift functions, which we consider 
as elderly notation in the powerful ADM approach. The potential in this 
work should satisfy only two essential conditions: matrix L is to be 
nondegenerate and the first class constraints are to fulfill the algebra of 
hypersurface deformations (necessary and sufficient conditions for this are 
equations (j57jl . fl58|) ). The special potential choice proposed by de Rham- 
Gabadadze,Tolley [I] will be treated in another article. 
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A First class constraints algebra in bigravity 

In two metric theory, as in GR, there are four first class constraints and 
their algebra is the algebra of hypersurface deformations. In computation 
of Poisson brackets between constraints TZ, IZi we treat u, u l not as canoni- 
cal variables, but as functions, because their conjugate momenta n u ,7i u i do 
not appear in the computation. Then potential U in these calculations has 
nonzero Poisson brackets with gravitational momenta , 7T 1 - 7 ' only, giving as 
results its derivatives with respect to induced metrics dU /drjij and dU /djij. 
First, consider equation fl55|) . generators %i and %i are mutually commuting, 
and each of them separately satisfies equation (13"T1) . therefore, equation (1531) 



is also fulfilled. L.h.s. of equation (|54l) can be given as follows 

{K{x),n(y)} = {n(x),n(y)}+ (71) 

+ {uH(x),uH(y)} + {u i H i (x),u j H j (y)}+ (72) 

+ {un{x),u i H i (y)} + {u i H i (x),uH(y)} + (73) 

+ {H(x),U(y)} + {U(x),H(y)} + (74) 

+ {uH(x),U(y)} + {U(x),uH(y)}+ (75) 

+ {u'n^^iy^ + iUix)^^)}. (76) 

Lines from 1 to 3 can be easily computed by using formulas (j30l) . fl3Tl) . fl32|) : 

= (v lk (x)n k (x) + v tk (y)n k ( y )) s ti (x, y ) + (77) 

+ u(x)u(y) (Y h (x)U k (x) + Y h (y)-H k (y)) 5 ti (x, y) + (78) 

+ u\x)u\y) (Hi(y)5, k (x, y) + U k {x)) 5 ti (x, y)) - (79) 

- u i (x)u(y)n(x)5, i (x,y) + u i (y)u(x)n(y)5, i (y, x), (80) 

line 4 and line 5 give zero results, as they are antisymmetric in x, y whereas 
each one contain 5-function. Last line gives a contribution 

dU dU 

Zu^jk- — (x)5 ti (x, y) - 2u 3 j jk - — (y)Si(y, x). (81) 

Taking into account relations 

f(y)fi,i( x > y) = f(x)S,i(x, y) + f ti 8(x, y), 5 ti (x, y) = -5,i(y, x), (82) 
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we can transform the result to the following form: 

{K{x),K{y)} = [Q\x) + Q\y)] S ti (x, y), (83) 

where 

dU 

= r) lk H k - uu l U + (Y V + u l u k )U k + 2u^ jk — , (84) 
To fulfill equation f )54|) it is necessary to have 

Q* =rj ik (H k + H k ). (85) 



As equation (J54J) is to be valid for Dirac brackets, we can use second class 
constraints in the expression obtained for Q\ i.e. replace %i by derivatives 
of the potential with respect for u, u\ so this requires to treat equation 
as a necessary condition for (|54p . 
At last, let us check equation (1561) : 

{Ki(x),K(y)} = {U i (x),H(y)} + 

+ {Hi(x),uU(y)} + {Ui^yUiiy)} + 

+ {n i (x),U(y)} + {ti l (x),U(y)} = 

= n(x)5 ti (x,y) + 

+ u(y)H(x)6 ii (x,y) + u 3 {y){Ui{y) + H j (x))5 ii (x,y) + 

/ . . <9?7 9C/ , A . , . 
+ 2 77 im (xJ- (x) +7j m (xJ- (x) 5 >n (x,y) + 



+ 



o 1 I till r\ I O ' 11111,1 r\ 



Proceeding in the same way to the previous calculation, and interpreting 
this bracket as Dirac bracket, we can use second class constraints in the 
expressions obtained. Then after exploiting formulas ( I8"2"j) . it is possible to 
demonstrate that in order to fulfill equation fl56|) the potential has to satisfy 
condition (l5j 



B Second class constraints and Dirac brack- 
ets in bigravity 



To make our presentation compact let us introduce notations: u a = (u,u l j, 
n a = (tTujTTuO) = CH>%) ? a — 1, -., 4. We denote eight second class 
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constraints as xa, — 1, ■•> 8, so now we have 

XA=\«a,na+j^y (92) 

Then matrix of Poisson brackets for these constraints has the following struc- 
ture 

where 

+ Q-^(x),H b (v) + (v) | , (95) 

we suppose that in the general case matrix L is invertable, then matrix 
||{/64)Xb}|| is invertable too, and its inverse is as follows 



h-\x)K.{x,y)L-\y) L" 1 (x)5(x,y) 
-L~ 1 (x)5(x,y) 



(96) 



Dirac brackets are given in the following way 

{F, G} D = {F, G}-JdxJ dy{F, X a(x)}C ab (x, y){ XB (y), G} (97) 

and in cases when both functionals F, G do not depend on variables u, u\ n u , ir u i, 
these brackets coincide with the Poisson ones. Therefore, if constraints are 
explicitely solved for u,u l , and these variables are replaced by constraint so- 
lutions, then for the rest of variables rjij, IP- 7 , 7^,71"^, ^4, IP 4 , 0a,tt a Dirac 
brackets coinside with canonical Poisson brackets: 

rrn ^ r j3 f5FSG 5F 5G 5F 5G 5F 5G _ \ 
{F,G} D = Jd*x{— — + — mx + — — + — _-(F«G)J. 

(98) 

But sometimes an explicit solving of constraints may be difficult, then one 
may write bigravity Hamiltonian by means of potential U, which is a function 
of variables u, u l and of the two induced metrics %-,7ij (ISBT) . In such a 
situation one has to use nontrivial Dirac brackets in obtaining Hamiltonian 
equations or in solving other problems: 

{j m n(x),u a (y)} D = -L- ia \y){ lmn (x),H b (y)}, (99) 
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{n™(x),u a (y)} D = L- ia \y) ( J*£—6(x,y) - {n mn (x), H b (y)}) , (100) 



du h d^n 



d 2 U 
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